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ABSTRACT 

In spite of the efforts made in the latest years, still there is no comprehensive ex- 
planation for the chemical anomalies of globular cluster stars. Among these, the most 
I striking is oxygen depletion, which reaches values down to [0/Fe]~ —0.4 in most clus- 

ters, but in M13 it goes down to less than [0/Fe]'^ —1. In this work we suggest that 
the anomalies are due to the super position of two different events: 1) primordial self- 
f~i , enrichment: this is asked to explain the oxygen depletion down to a minimum value 

Oh' [O/Fe]'^ —0.4; 2) extra mixing in a fraction of the stars already born with anomalous 

composition: these objects, starting with already low [0/Fe], will reduce the oxygen 
^ . abundance down to the most extreme values. Contrary to other models that invoke 

^ ' extra mixing to explain the chemical anomalies, we suggest that it is active only if 

^ I there is a fraction of the stars in which the primordial composition is not only oxy- 

gen depleted, but also extremely helium rich (Y~ 0.4), as found in a few GCs from 
their main sequence multiplicity. We propose that the rotational evolution (and an 
associated extra mixing) of extremely helium rich stars may be affected by the fact 
that they develop a very small or non existent molecular weight barrier during the 
I evolution. We show that extra mixing in these stars, having initial chemistry that has 

, already been CNO processed, affects mainly the oxygen abundance, and to a much 

ff^ ■ smaller extent if affects the sodium abundance. The model also predicts a large fluo- 

rine depletion concomitant with the oxygen depletion, and a further enhancement of 
the surface helium abundance, which reaches values close to Y=0.5 in the computed 
' models. We stress that, in this tentative explanation, those stars that are primordially 

, 0-depleted, but are not extremely helium rich do not suffer deep extra mixing. 
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1 INTRODUCTION (Hot Bottom Burning -HBB) in low velocity winds, that 

, . r , , , ^, , may remain into the cluster and form the second generation 

the recent observations 01 abundance spreads among Glob- i„ — rn — — 7; — . 11 ITT — I 1 — nUr^mh 

% /rrr^, (^Cottrell fc Da Costall 19811 : IVentura et al.ll2001l l. 

ular Clusters stars, now observed also at the turnott (iO) . , r r , , ,■ • , 

. , , — ; — n u^^j v , u) Wl^KMb: wmds ot last rotatmg massive stars, ihe 

and among the subgiants (e.g.. IGratton et al.l 120011 ) show . ,. , , ■, r n , 

that these anomalies must be mainly attributed to some 

process of "self-enrichment" occurring at the first stages of 



primordial gas cloud, from which the first generation stars 
are born, is contaminated with the winds ejected from lastly 
^ . , rotating massive stars during the first few million years 

the life of the cluster. Two mam sources have been proposed ^^^^^^^ jMaeder fc Mevnedbood: lOecressin et al.l 

as basis tor the selt-enrichment: U^^A rr^, ■ , , ■ 7- , , r , 

120061). Ihis model implies a very last second stage oi star 

(i) WMAGB: Winds of massive AGBs. During the epoch formation, preceding the supernova events whose energetic 
starting some ~ 5 x lO'^yr from the birth of the first stellar input on the cluster gas may preclude any further star for- 
generation, when all the Supernovae have already exploded mation. 
(carrying away from the clusters their high velocity ejecta) 

the massive Asymptotic Giant Branch (AGB) stars evolve, In the case of WFRMS, the primordial gas from which the 

losing their envelopes cycled through the hot CNO-cycle fist stellar generation is born is certainly still present when 

the nuclearly processed stellar envelopes are ejected, and 

these latter may be diluted at different degrees in the primor- 

^ „ J i ■ r i r-wT^\ i dial gas before new stars are formed. In the case of WMAGB, 

E-mail: dantonaiaoa-roma.mat.it (FU); ventura@oa- ° ' 

roma inaf it (PV) remnant gas of the first generation may still be present and 

t This work has been supported through PRIN INAF 2005 "Ex- dilute the ejecta either if the polluting stars are only the 

perimenting stellar nucleosynthesis in clean environments" most massive AGBs (those referred to as super- AGBs) , or 
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if, e.g., GCs are formed into dw arf galaxies at the time of 
the evolution of massive AGBs l|Bekki et al.ll2007l ). In the 
latter case, the gas ejected from the AGBs may be collected 
and pollute the gas present in the region where the cluster 
is being formed. Otherwhise, it is scarcely credible that the 
winds are diluted with residual gas before forming stars of 
the second generation (jCottrell fc Da CostairToSlI ). as sev- 
eral decades of million years have passes, and the supernovae 
have cleaned up the cluster. In fact, it has also been sug- 
gested that the spreads in chemical abundances are actually 
due to the birth of successive generation of stars directly 
from the ejecta of the massive AGBs of the first generation 
l|D'Antona et al.ll2002l '). 

Neither model has been yet explored from a dynamic 
point of view. The role of the fast post-main sequence 
winds and of their energetic input in the case of WFRMS 
is not clear. An escape from the fact that the supernova 
events will be dramatic for the gas environment in the clus- 
ters is to consider that stars more massive than ~25 Mq 
may end their life collapsing into black holes. The advan- 
tage of this scenario is that the requirements on the initial 
mass budget in massive star s is less severe than the con- 
straints of the AGB scenario llPrantzos fc Charbonnelll2006l : 
iDecressin et aili200(t ). In regard to the chemistry, both mod- 
els are still not fully adequate, although there is scarce doubt 
that the abundance anomalies we see are born in matter pro- 
cessed through the hot CNO cycle, wit h little -if any- influ- 
ence of the third dredge up in AGBs l| Ventura et ahllioo^ : 
IVentura fc D'Antona|[2005al '). 

An important prediction of both scenarios is that the 
second stellar generation must have helium content en- 
hanced, at some degree, with respect to the helium con- 
tent of the first stellar generation, that probably refiects the 
Big Bang abundance. In recent years, the prediction of the 
spread that is naturally associated with the helium abun- 
dance in the AGB ejecta, allowed to interpret the morphol- 
ogy of extended HBs in some globular clusters in terms of 
a spread in the initial helium content of the cluster stars 
llD'Antona et al.l |2002| : iD'Antona fc Caloil '2004'; 'Lee et al' 



Interestingly enough, after the ID'Antona fc Caloi 



interpretation of the HB of NGC 2808 in terms of two 
populations differing in initial helium content, a very pecu- 
liar helium distribution in the cluster stars was inferred from 
the presence of a "blue main sequ ence", containing ^15 - 
20% of stars, in the same cluster (|D'Antona et ahlbOOSl) , 
similar to w h at had been found in uCen l|Bedin et al. 1 12004 
Norrid |2004|: IPiotto et all 120051 ). New and better data by 
all 



Piotto et al.l (|2007l ) have actually shown that the main se- 



quence of NGC 2808 actu ally is formed by three main com- 
ponents, precisely as the ID'Antona et al.l (|2005h interpre- 
tation of the HB population required (see their Figure 9). 
In general, the main sequence morphology can allow to de- 
tect differences in heliu m content only when the helium con- 
tent is really very large (|D'Antona et al.|[2002l : ISalaris et al.l 
l2006f ). In fact, both in NGC 2808 and tjCen, a straightfor- 
ward interpretation through stellar models assigns Y~0.4 
to the bluest main sequence. This population however is 
only the "tip of the iceberg", as the interpretation of the 
HB morphologies requires in many clusters that 50% or 
more of the stars have enhanced helium abu ndance (e.g. 
ICaloi fc D'Antona|[2007l : ID'Antona et all 120061 ). The most 
straightforward interpretation of existing data is that most 



clusters have a population with enhanced helium (up to 
Y~0.3-0.33), but only a few have a small population with 
very htgh helium, such as NGC 2808. We now try to under- 
stand whether we can recognize a similar behaviour also in 
the sodium versus oxygen anticorrelation, which re presents 
the m ost clear evidence of chemical anomalies (e.g. ISnedenI 
Il99ll ). This anticorrelation is present in most clusters ex- 
amined so far. The maximum reduction found for oxygen 

in most clusters reaches a minimum value [O/Fe] 0.43, 

while the "normal" popul ation in GCs and the halo stars 
have [O/Fe]-- 0.3 (see, e.g. ICarretta et al.|[2006l . their Fig. 
5). Sodium in the stars with reduced oxygen is in the range 
[Na/Fe]~ 0.3 - 0.7, while [Na/Fe]-0 in the "normal" stars. 
The Na-O anticorrelation is al so present among th e TO 
and subgiant stars, as shown bv lGratton et al.l l|200lf ). and 
therefore it can not inputed to 'in situ' mixing. In addi- 
tion to these typical anomalies, some clusters have a small 
number of stars showing more extreme behaviour: while the 
largest sodium abundance remains at ^^-1-0.7, oxygen can be 
reduced d own to [0/Fe]~ — 1.2 in a few stars in the cluster 
M13 (e.g. iKraft et al.lll993l ). The cluster NGC 2808, which 
shows the "blue main sequence" with Y~0.40 comprising 
about 15% of the cluster stars, also seems to have a sim- 
ila r fraction of ex t remel y oxygen depleted stars, according 
to ICarretta et al.l (|2006l). We are led to suggest, purely on 
observational grounds, that the stars having the lowest oxy- 
gen values belong to the same group of those stars having 
the highest initial helium abundance. Further, until today, 
the most extreme anomalies are only found in red giants. 
In this work, we try to understand whether these informa- 
tion can be used to derive some simple hint on the inter- 
pretation of the abundance anomalies, to be used as a con- 
straint on the formation of the second stellar generation. 
In Section 2 we examine the possible schemes for the ex- 
planation of the Na-O anticorrelation. Among the different 
possibilities, we suggest that some among the "typically ex- 
treme" ([O/Fe]'-- 0.4, [Na/Fe]— 1-0.5-0.7) stars evolve into 
the most extreme Oxygen poor giants due to an "in situ" 
extra mixing. In other words, we will make the case that 
not all the abundance anomalies can be attributed to the 
chemistry processed in the progenitor stars (being either 
WMAGB or WFRMS poUutors), and that we need some 
process of extra-mixing in the red giant stage to explain the 
most extreme anomalies. In some clusters only, we suggest 
that there is extra mixing in red giants, and correlate it to 
the existence of an extremely helium rich population. We 
justify this suggestion on the basis of the different model 
structures during the RGB stage for differing stellar helium 
content (in Appendix). In Section 3 we summarize the con- 
cept of "canonical" and "non canonical" extra-mixing in 
the literature. In Section 4 we model parametrically extra- 
mixing during the RGB evolution. We show that the results 
of extra mixing depend on the abundances assumed for the 
initial stellar chemistry. Starting from already anomalous 
composition, deep extra mixing will mainly affect the oxy- 
gen abundance. Table 1 summarizes the possible scenarios 
to understand the GC chemical anomalies. We also show the 
effect of this parmetric extra mixing on the Fluorine abun- 



1 as usual [X/Fe]=log(Nx/Nj7e) - log(Nx/NFe)Q , where N in- 
dicates the abundances in number. 
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Figure 1. Left: observed data that define the anticorrelation sodium vs. oxygen in the stars of several GCs. stars: M13, open squares: 
M3 (both from Sneden et al. 2004); full squares: NGC 6752 from Grundahl et al. 2002; full triangles: M4 and open triangles: M5 (both 
from Ivans et al. 1999); full dots: NGC 2808 from Carretta et al. 2006. The three curves represent mixtures of different quantities of the 
two compositions represented by the squares at each extreme. Going from left to right, the composition at the extreme left is diluted 
by 0, 10, 20 ... 90% of mass having the composition at the extreme right. Right: this figure shows only the M13 (stars) and MS (open 
squares) data, and the possible dilution curve for M13. The (red) dash-dotted line represent the result of extra mixing from the sequence 
4 of Table 1. 



dance, and point out that the models, starting already with 
Y=0.4, show a surface helium increase up to Y~0.5, which 
could be important to understand the very hot subdwarfs 
present in the HB of some clusters. 



2 EXPLAINING THE OBSERVATIONAL 
SODIUM VS. OXYGEN 
ANTICORRELATION 

2.1 The abundance patterns and the dilution 
model 

We show in Figure [T] the Sodium versus Oxygen abun- 
dances among the stars in several GCs. In contrast with 
the halo field stars, which occupy only the lower right side 
of the figure (high oxygen, normal sodium), m any of the GC 
stars show low oxygen and higher sodium (e.g. ISneden et al.l 
[2004). The case of the clusters M3 and M13 is shown sep- 
arately in the right side of the figure. The continuous lines 
with squares show a possible and straightforward interpreta- 
tion of the abundance patterns: if we mix gas having the high 
Na-low O composition of the square at the extreme left with 
"normal" matter having the low Na-high O composition rep- 
resented by the square at the extreme right, we obtain all 
the intermediate compositions. The matter giving origin to 
the second stellar generation should then have [0/Fe]~-1.2 
and [Na/Fe]~0.7 in the case of M13, or less extreme values 



for other clusters, e.g. [O/Fe]~-0.4 and [Na/Fe]~0.3 for the 
case of M5. 

In addition to M13, extremely small oxygen abundances 
are possibly present also am ong the stars in NGC 2808, as 
shown by the upper limits of lCarretta et al] (|200Q ) analysis. 
For M13 and possibly for NGC 2808, we seem to need a 
site of nucleosynthesis (the envelopes of massive stars, or of 
AGB stars) in which the oxygen abundance is reduced by 
~1.5dex. 



2.2 Models of FRMS 

The models bv lDecressin et ahl (|2006l ) show that the matter 
lost in the winds has composition which indeed keeps the 
signatures of CNO processing. For what concerns the Oxy- 
gen abundance, however, the maximun reduction they find 
is from an initial abundance of 3.0x10"'* to 1.1 x 10~* for 
a 120 Mq evolution (their Table 5). At the same time, the 
largest oxygen anomalies are found together with extreme 
helium overabundances (Y up to 0.5 in the quoted 120 Mq 
evolution). These overabundance are of great importance, 
as we will see later on. Thus, starting from an a-enhanced 
[O/Fe]~0.3 for the oxygen of the first stellar generation, the 
minimum abundance which could b e found in second gen - 
eration stars would be [O/Fe]~-0.4. iDecressin et all (|2006l ) 
also consider the fact that the Lithium abundances so far ob- 
served for turnoff stars in GCs are smaller than the primor- 
dial Lithium abundance, but still quite large. They use the 
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Table 1. Possible scenarios: composition of the polluting gas 



Source 


Cluster 


[O] 


[Na] 


Y 


mechanism 


? 


WFRMS 


extreme GCs: 


-1.3 


0.7 


0.40 


dilution with 1st gen. gas 


not reliable: pollution by 




(M13 ; NGC 2808) 










fast winds is necessary 


WFRMS 


extreme GCs 


-0.4 


0.5 


0.40 


dilution with 1st gen. gas, 














extramixing for extreme c.c. 


proposed here 


WFRMS 


less extreme GCs 


-0.4 


0.5 


<0.35 


dilution — no extra mixing 


proposed here 


WMAGB 


extreme GCs 


-0.4 


0.5 


0.40 


undiluted gas, or dilution with 1st gen. gas 


proposed here 



extramixing for extreme c.c. 
WMAGB less extreme GCs -0.4 - 0.0 0.3 - 0.5 <0.35 undiluted gas - no extra mixing 



value AfLi') = 1.9j^bv |Pasquini et all (120051 ) to infer that the 
matter processed in the FRMS must be diluted with pristine 
matter (having A(Li) =2.61, Coc et al. 2004) in the propor- 
tion 2/3 from WFRMS and 1/3 from the pristine matter, 
in order to preserve this relatively low Lithium abundance. 
In this case, the minimum oxygen value becomes [O/Fe]^- 
0.3. However the data on Lithium are scarce enough that 
we should postpone any judgment of this model based on 
Lithium to a more complete observational situation. There 
is an alternative possibility to reduce the oxygen content, 
which is worth mentioning: if we do not consider the average 
values for the abundances in the winds of FRMS, but their 
time evolution, much smaller oxygen abundances (down to 
~ 7 X 10~®) can be achieved (Figure 6 in Decressin et al. 
2006). If the winds may "directly" form stars, some of them 
might be as extreme as required by M13 observations. No- 
tice however that, in this case, we need to use the chemistry 
which the winds reach in the phase of He-burning. At this 
stage, the stars are no longer fast rotating and the winds 
have high velocities, so that it is difficult to understand how 
they can be kept into the GC. 



2.3 Models for massive AGB stars 

For certain oxygen abundances as low as [0/Fe]~-l are not 
achieved in the yields of intermediate mas s AGBs compute d 
until today ([Vcntura & D'Antona 2006; Bekki et al.ll2007l ). 
Nevertheless, recent modelling of intermediate mass AGB 
stars of low metallicity is not doing a bad job to provide 
a reasonable framework to explain a gre at part of che GC 
chemical anomalies. Ventura & D'Antona'2005al: l2005bl have 
considered carefully the role of convection model and nuclear 
reaction rates in fulfilling the requirements of the chemical 
patterns of abundances to be produced. In particular, they 
showed that the Na-O anticorrelation requires very efflcient 
convection, to process the CNO through the ON cycle by 
HBB. A problem arises as, at the same temperature nec- 
essar y for destroying oxygen, also Na is efliciently burned 
fe.g. iDenissenkov fc Herwigj |2003| ). After some very para- 
metric attempts to solve this problem by means of a very 
small overshooting associate d to the thermal pulse phase 
l|Ventura fc D'Antona|[2005d ). a possible, very stable, solu- 
tion is obtained by lowering the reaction rate ■^''Na(p,a)^''Ne 
within limits acceptable in the allowed limits for this cross 



section , and raising by a factor ^ 2 the initial neon abun- 
dance l| Ventura fc P'Antonal [200^ ). Nevertheless, the most 
recent finding is that, for M> 5Mq, there is a sort of "sat- 
uration" in the combined effect of HBB and mass loss, and 
the oxygen reduction can not go beyond ~ — 0.7dejjf| (Ven- 
tura and D'Antona 2007, in preparation). The sodium en- 
hancement does not go beyond ~ -|-0.5dex, but this may be 
linked to the initial values assumed for the abundance of 
^■^Ne. These models are then able to explain the bulk of the 
Na-0 anticorrelation, which is present in most GCs, but not 
the "extreme" stars of M13 and NGC 2808. 



2.4 A closer look at the data 

From 2.2 and 2.3 we conclude that both the possible sites 
for the nucleosynthesis necessary to explain the chemical 
anomalies are not consistent with the very low oxygen abun- 
dances of some giants of M13. Let us examine these in detail: 
the data provided bv ISneden et al.l (|2004l ) for M13 contain 
8 stars with [O/Fe]<-0.4 and [Na/Fe] from ~0.4 to ~0.8; 
among these, 6 are brighter than Mv=-2.3 and are there- 
fore essentially at the tip of the giant branch. Beside, 5 
of them are surely vari able, either semi- or irregular (1-48, 
11-90, L835, L940 from lKopacki et al.| l2003'. L954 in JWeltvl 
il985); one is a suspected variable (11-67 in .Wclty , l,98j). For 
the two lower luminosity stars (III-59 and 1-12) there are no 
indications of variability. The extreme abundances and (very 
low) upper limits for oxygen con cern a fraction of the stars in 
NGC 2808 (|Carretta et a l. 20061. b ut also these data refer to 
evolved giants. The compilation bv lCarretta et al.l l|2006l ) of 
the available data for the Na-O anticorrelation shows clearly 
that the scarcely evolved stars in GCs do not show abun- 
dances smaller than [O/Fe]~-0.4 and that the data below 
[O/Fe]=-0.4 all belong either to M13 or to NGC 2808. In 
addition, although u! Gen is certainly much more compli- 
cated than a typica l GC, the sample of giants examined by 
iNorris fc Da Costal (|l995h contains a few upper limits on 
oxygen abundances below [0/Fe]~ —0.4, for CN weak stars 
with high sodium abundances. From this discussion, it ap- 
pears then possible that in some clusters, for a fraction of 
the evolved stars, an in situ mechanism acts to produce the 
extreme abundances. 



2 the notation is A(Li)=log [N(Li)/N(H)]+12 



^ starting from [0/Fo]~ 0.3, this means that in these models 
[O/Fe] > -0.4, as in the case of WFRMS. 
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2.5 Are the stars with extreme oxygen depletion 
also very helium rich? 

Let us assume that the very helium rich sequences are 
present only in the clusters having the most extreme anoma- 
lies. Then the minimum helium co ntent of the processe d 
matter must be Y~0.4. As shown bv lDecressin et al.l(|2006h . 
the winds of massive rotating stars are heavily enriched in 
helium. In the dilution model, therefore, the most extreme 
anomalies — those associated with the pure stellar ejecta — 
would also correspond to the stars having the most extreme 
helium abundances. For the WMAGB model, the evolution 
of AGBs is characterized on the one side by HBB, and 
on the other side by an efflcient second dredge up of he- 
lium, which raises the helium abundance in the envolopes 
to values from Y~0.27 (for a 4M0Star) to Y~0.3-0.32 for 
a 6Mq. There are until now no models computed for the 
super-AGBs (masses 7-8 M©) of population II, and we can 
hypothize that these would reach the extreme helium re- 
quired, Y~0.4. Therefore, both models in principle suggest 
that the maximum helium values are linked to the most ex- 
treme anomalies. From an observati onal point of vi e w, let 
us remember an intriguing result bv ICarretta et al] (|2006l ) 
concerning NGC 2808. Dividing the cluster giants into three 
groups, with decreasing [0/Na], they notice that there is 
a very small increase in [Fe/H] in the three groups, ex- 
actly as one would expect if the Na-O anomaly is associ- 
ated with a larger and larger helium content of the studied 
atmospheres: thus it is possible that the extremely low oxy- 
gen giants in NGC 2808 belong to the same group of stars 
which populate the blue main sequence, and the extreme 
blue tails of the horizontal branch distribution. Unfortu- 
nately, in M13, the only other cluster showing extremely 
low oxygen abundances, an high precision photometry of the 
MS is lacking, to see whether a blue main sequence can be 
isolated. For a; Cen, the other cluster showing a blue main 
sequenc e which can be interpreted in terms of very hig h 
helium (iBedin et al.ll2004l : lNorrisll2004 iPiotto et al]|2005l ). 
there might also be giants with extreme valu es for the an- 
ticorrelation Na-O (|Norris fc Da Costal 1 19951 ). but further 
studies are needed. 

We then have the following problem: the models may 
probably, in some cases, provide ejecta with extremely high 
helium content, but section 2.2 and 2.3 suggest that this 
same matter can not have the extreme oxygen content of 
M13 giants. These very low abundances, then, should be 
due to some in situ mechanism (extra mixing) taking place 
only in these stars. An important question we should ask 
is whether there is any physical reason to link a very high 
initial helium abundance to extra-mixing during the RGB 
evolution. We give a speculative answer in Appendix A. 



3 CANONICAL AND NON CANONICAL 
EXTRA-MIXING IN RED GIANTS 

Standard mixing due to the convection description based 
on the Schwarzschild instability criterion is not able to 
explain the decline of the abundance of Carbon and the 
strong reduction of the ^^C/^^C isotopic ratio in upper 
red giant br anch (RGB) stars, both in the field and in 
star clusters (|Snedenlll99lh . and the decrease in surface Li 



abundance with increase d luminosity l|Pilachowski et akl 
1 19931 : iGratton et aLllioOol ). Additional mixing seems to be 
efficient in the evolving stars, as soon as the H-burning 
shell, advancing in mass, erases the chemical composition 
discontinuity left behind by the bottom of convective 
envelope at the end of the first dredge-up (|Gratton et al.l 
I2OOOI : IShetronell2003l '). At this moment, the evolution of red 
giants slows down for a while (see Fig. IA3|) . This produces 
a prominent feature (the RGB bump) i n the differential lu- 
minosity function s of globular clusters (|Zoccali et al.|[l999l : 
iRiello et al.ll2003l ). The presence of the RGB bump can be 
taken as an experimental indication that the "standard" 
models of red giants actually make a good job in describing 
the evolution below the bump luminosity, and that no 
mechanisms of extra mixing in low-mass stars are acting, 
on the lower RGB, otherwhise the discontinuity left by 
convection wo uld be smeared out an d the bump would not 
be produced l|Chaname et al.l l2005t ). The most credited 
hypothesis for preventing extra mixing below the bump 
is that the H-burning shell is shielded against mixing 
by the gradient in the mean molecula r we ight associated 
with the c omposition discontinuitv JSweiga rt fc Mcngdl 



19791: 



Charbonnel. Brown fc WallersteinI 



Denissenkov fc VanderiBerg 20031) or, in 



opera tes very slowly (IChaname et al. 
'2OO6!). Very recently. lEggleton et al 



I2OO5I: 



any case. 



Palacios et al.l 



(|2006l ) proposed a 



different analysis: mixing arises between the supposedly 
stable and radiative zone between the hydrogen-burning 
shell and the base of the convective envelope, due to 
Rayleigh- Taylor instability within a zone just above the 
hydrogen-burning shell, where the burning of the "^He left 
behind by the retreating convective envelope lowers slightly 
the mean molecular weight, leading to a local minimum. 

For the majority of upper RGB stars, the surface abun- 
dance anomalies require that this in situ mixing penetrates 
down to the outer part of the H-shell, where the CN branch 
of the CNO-cycle is operating . This is the "canonical ex- 
tra m ixing" , in the definition of iDenissenkov fc VandenBergj 
'2003!. These authors, parametrizing the depth and rate (dif- 
fusion coefficient) of extra-mixing below the formally convec- 
tive region, find that the values which describe the "canon- 
ical" extra mixing are within the range from A logT ~ 
0.19 and Dmix~ 4 x 10* cm^s"^ to Alog T ~ 0.22 and 
Dmix~ 8x 10* cm'^s"^. Here A logT is the difference between 
the logarithms of temperature at the base of the H-burning 
shell and at the maximum depth of extra mixing. These val- 
ues then can explain the evolutionary decline of [C/Fe] in the 
upper RGB stars, both in the field and in star clusters. How- 
ever, canonical extra mixing does not affect O, Na, Mg and 
Al, the elements which show anomalous abundances in the 
stars of GCs. A "non canonical" extra mixing in upper RGB 
stars, penetrating the H-burning shell deeper than in the 
canonical case, c ould dredg e up mat erial deficient i n O and 
enriched in Na (|Denisenkov fc Denisenkoval Il990l ). These 
models constituted an alternative view to the problem of GC 
chemical anomalies, until it was clear that these anomalies 
were present also in scarcely evolved stars. Once established 
that the self-enrichment due to star formation in matter 
either from WFRMS and/or from WMAGB is a better sce- 
nario to explain globally the chemical anomalies in GC stars, 
it remains to be understood whether mixing plays any ad- 
ditional role. The possibility of enhanced, faster and deeper 
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Figure 2. The modification of ^ONe, ^^jsjj, and ^^Na profiles inside a red giant having solar scaled initial ^^Ne (figures on the left) or 
depleted ^^Ne (figures on the right). In the top panels, the upper curves (almost unchanged) refer to the ^^Ne abundances, and the lower 
ones to the ^^jsje. Alternate dotted and full lines are adopted to indicate successive evolutionary phases. The fractional mass is given as 
abscissa, the total mass is O.64M0. In early phases of the RGB evolution, the high ^^Ne models burn it to ^^Na producing a profile of 
high sodium inside the giant. If the ^^Ne is already low, the ^^Na does not increase. 



Table 2. Evolution of chemical abundances in models with Di,iix= 5 X lO^'cm^s ^ 
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0.05 


4 f 
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0.06 


5 ml 


600 


0.06 


6 f 
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0.06 





["N] 




-0.862 


0.904 


-0.208 


-2.009 


0.936 


-0.209 


-1.353 


1.128 


-1.404 


-1.373 


1.119 


-1.140 


-1.356 


1.129 


-1.467 


-1.393 


1.096 


-0.811 



[19F] 


[22 Ne] 


[23Na] 


-1.104 


-1.367 


0.539 


-1.106 


-1.598 


0.548 


-2.978 


-2.336 


0.761 


-2.740 


-2.388 


0.687 


-3.047 


-2.358 


0.840 


-2.378 


-2.416 


0.672 



mixing, has not been neglected as an add itional mechanism 
during the RGB evolution: e.g. ;Denisscnkov &: VandenBerd 
l|2003l ) have proposed that in some upper RGB stars canon- 
ical extra mixing may become deeper, driven by the dif- 
ferential rotation of the stellar radiative zones due, e.g., 
to spinning up of clo se binary members as a result of 



tidal synchronization ( Denissenkov et all l2006al ). In addi- 
tion, iDenissenkov et al.l l|2006bD proposed enhanced extra 



mixing to explain the Fluorine abundance variations in the 
red giants of the cluster M4. 

In the following we will assume that indeed a non 
canonical extra mixing is ejjicient only in most helium rich 
fraction of cluster stars, and we model its effect on the 
surface abundances. This combination o f a primordial and 
deep m ixing scenario has been invoked bv lDenissenkov et al.l 
l|l99d ) to explain the magnesium vs. aluminum anticorrela- 
tion in GC red giants, and we follow their line of reasoning, 
although we will devote to another study the comparison of 
the magnesium and aluminum abundances, and discuss here 



only the sodium and oxygen abundance behaviour. There- 
fore, we will assume that the initial abundances of the evolv- 
ing giants are consistent with the chemistry obtained either 
in the WFRMS or in the WMAGBs. In order to have a 
precise choice for the computation of the models described 
below, we adopt the chemistry of the e.jecta of a 5 sta r 
according to our AGB models (jVentura fc D'Antonair2006l '). 
but we assume the high helium abundance which would re- 
sult from a higher mass (Y=0.4). 

3.1 Models for non canonical extramixing 

The models are computed by making use of the ATO N code , 
described in its latest version in Ventura & D ' Antona l2005al : 
[2005b. We compute standard models, in which mixing is 
purely due to convection, and no kind of overshooting is al- 
lowed below the external convective region. A description 
of these models is given in Appendix A. Here we describe 
the models with extra mixing, which refer to M=0.64 Mq 
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Table 3. Evolution with different Dmixand starting phase 
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Figure 3. Sodium (top) and Oxygen (bottom) evolution in more 
evolved models of the 0.65 Mq. The dotted lines refer to stan- 
dard models, the continuous lines refer to models with deep extra 
mixing. Oxygen becomes more and more depleted at the shell 
boundary, by conversion to ^^N, so that fast extra mixing is at 
once effective in reducing its surface abundance. On the contrary, 
only at high luminosity further burning of 22 Ne finally establishes 
a profile of 23Na into the star. Extra mixing can raise the surface 
abundance of 23Na only in very late phases. 



with initial helium mass fraction Y=0.4, and metallicity in 
mass fraction Z=0.002. The helium abundance is selected 
following our reasoning which links the extra mixing to the 
initially very high helium abundance both as a consequence 
of observations (Section 2) and of the interpretation of the 
RGB evolution as a function of Y (Appendix A). As we will 
justify in the Appendix, we do not link the occurrence of ex- 
tra mixing to a particular stage of evolution along the RGB, 
but we consider different initial stages, from luminosities 
at the beginning of the RGB evolution, up to luminosities 
larger than those of the (very mild) RGB bump for the cho- 
sen Y. The mass 0.64 Mq is chosen as representative of a 
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log L/L< 



Figure 4. Oxygen and sodium evolution as a function of the 
luminosity in the models of Table 2. 



mass presently evolving in a typical GC (age of ~ 12Gyr) 
for Y=0.40 and Z=0.002. 

Some model results are listed in Table 2 and 3. At this 
stage, we do not aim to understand what is the specific mech- 
anism for extra mixing, but only which are the consequences 
for the surface abundances of the giants presently evolving 
in GCs. We adopt the description by Denissenkov & Van- 
deuBerg (2003), and assume those depths and rates of extra- 
mixing, below the formally convective region, necessary to 
obtain sensible changes of the surface abundance of oxygen. 
Values of AlogT ~ 0.05 - 0.06 and Dmix~ 0.1 - 5 x lO'' 
cm2s~^ are adopted. We also followed an evolution by con- 
sidering the mass loss by stellar wind. Table 2 lists in the first 
column an identification sequential number and the phase 
of evolution is represented. Sequence 1 ('st') is the sequence 
computed without extra mixing, while 'f ' means that we pro- 
vide the abundances at the end of the computed sequence, 
corresponding to the onset of the helium fiash. In the 'ml' 
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Table 4. HR diagram location and helium core of models along 
the standard 0.64 Mq sequence 



N 


logL/Lg 


log T,ff 


M(H)/Mq 


500 


0.567 


3.752 


0.099 


600 


0.909 


3.727 


0.165 


900 


1.604 


3.692 


0.234 


1050 


1.797 


3.681 


0.253 


1500 


2.147 


3.660 


0.291 


2000 


2.427 


3.641 


0.319 



sequence 5 of Table 2, mass loss is included in the compu- 
tation, according to the description by Reimers': 

M = 4.83 X lQ~^r]RLR/M (1) 

where L, R and M are respectively the stellar luminosity, 
radius and mass in solar units and rm is the free parameter. 
In this case we take 77jj=0.29. 

The third column lists the initial model at which extra 
mixing starts, beginning from the sequence with no extra 
mixing. The luminosity, Tea and helium core mass corre- 
sponding to the model numbers are listed in Table 4. The 
following columns provide the logarithm of abundance of 
the listed elements, with respect to the standard solar-scaled 
abundances. The values at model 1 of the standard sequence 
1 do not change along the evolution, so these element abun- 
dances also represent the starting abundances of the other 
sequences. 



4 CHEMICAL STRATIFICATION INSIDE THE 
RGS AND THE EFFECT OF THE INITIAL 
CHEMISTRY 

The effect of deep mixing depends on the profiles estab- 
lished within the star for the elements into consideration. 
These profiles depend on the one side on the temperatures 
reached within the inner layers during the evolution, and on 
the other side on the initial abundances of these same ele- 
ments and of those which produce them. For the elements 
which are our main concern now, in particular, oxygen and 
sodium, in fact, we wish to study the effect of deep ex- 
tra mixing starting from the initial mixture that comes out 
from the evolution of a 5 Mq of metal abundance in mass 
fraction Z=2xl0~'^, a ccord ing to the models described in 
IVentura" fc D'Antonal l|2006l ). These abundances are listed 
in the first line of Tables 2 and 3. In particular. Oxygen 
is already depleted by -0.508 dex with respect to its initial 
abundance ([O/Fe]=0.3fl, Sodium is already enhanced by 
+0.54dex, and ^'^Ne is depleted by -1.37dex. 

These initial abundances have an obvious but subtle 
effect on the result of extra mixing. In particular, let us 

* The choice of an initial abundance [0/Fe]~ —0.2 may be con- 
fusing, as we have stated that abundances as low as [0/Fe]~ —0.4 
may be achieved in massive AGB models. We preferred to use 
"directly" the abundances coming out from a precise stellar evo- 
lution, namely that of the 5Mq, as we have not yet completed 
the study of the evolution of more massive models which have 
the maximum oxygen depletion (Ventura and D'Antona 2007, in 
preparation) . 



first consider the Ne and Na abundances. Let us con- 
sider models starting from solar scaled abundances: The 
■^■^Ne abundance is initially almost 4 times larger than the 
■^■^Na. There is ample space for ^^Ne to be burned inside 
the giant and establish a large gradient of '^'^Na across the 
burning shell, even at relatively low luminosities. We show 
in the left side of Figure [2] the formation of this sodium 
gradient at several stages during the standard giant evolu- 
tion of a 0.64 Mq having solar scaled initial abundances of 
these two elements. The advancing shell burning shifts to 
larger fractionary masses the formation of high sodium and 
the burning of ^'^Ne by proton captures. But in the AGB 
envelopes which we assume to be progenitors of the giants 
which suffer extra mixing, already a great fraction of ^^Ne 
had been converted to ^'^Na by HBB. Therefore, during the 
first phases of RGB evolution, no further ^^Na is produced, 
and no sodium profile is established in the star, contrary 
to the case of abundances in solar ratios (right side of Fig- 
ure [2]). It is only later on during the RGB evolution, when 
the shell temperature is higher, that ^'^Na is again produced 
at the expenses of the burning of a very small fraction of 
the abundant ^"Ne. This is shown in the top of Figure O 
sodium enhancement at the border of the burning shell can 
be obtained only when the core mass becomes larger than 
~O.3M0, that is at high luminosity. In this figure the dot- 
ted lines represent the sodium (top) and oxygen (bottom) 
profiles in the standard model evolution. The abscissa, as 
usual, is the fractionary mass. The location of the abun- 
dance gradients corresponds to the core mass (and therefore 
luminosity -see Table 4) achieved at each model. We see 
that, while sodium is produced and oxygen depleted inside 
the star, nothing happens to the surface abundances if there 
is no extra mixing. The continuous lines show the evolution 
of the chemical profiles into some models of sequence 4 of 
Table 2. We see that oxygen begins to be depleted at much 
smaller core masses (luminosities) than those at which some 
effect on the increase of sodium begin to occur. Notice that 
we have started from an already depleted oxygen abundance: 
if we had started from solar scaled -or, better for standard 
population II, enhanced- oxygen, the dilution effect would 
have been much larger. In fact, in the classic interpretation 
of abundance anomalies by deep mixing, starting from stan- 
dard initial abundances is possible to interpret the whole 
sodium v s. oxygen anticorr elation in term of deep mixing 
(see, e.g.. IWeiss et al.ll2000l . their figure 9). 

The fact that the effect of extramixing is more impor- 
tant for the oxygen depletion than for the sodium increase 
is very clear by looking at Figure 2] which shows the abun- 
dances of oxygen and sodium vs. the stellar luminosity, for 
the extra mixing models detailed in Table 2. One of these 
results (model 4) is shown in the right side of Figure [T] 
We see then that in fact the very low oxygen abundance of 
M13 giants can be consistent with an in situ extra mixing 
starting at non extreme abundance anomalies ([O/Fe]^-0.4, 
[Na/Fe]'--^+0.5) provided by the massive AGB models. 

4.1 The Fluorine abundance 

Fluorine is already depleted by a factor ~10 by HBB in the 
5 Mq evolution from which we take the initial abundances. 
This depletion is consistent with the data by ISmith et al.l 
l|2005l ) for the giants of M4. This cluster should not be as 
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Figure 5. The modification of ^^F abundance for various evolutions (Table 3) are shown as a function of the oxygen abundance (left) 
and or the stellar luminosity (right figure). 



Table 5. Final helium mass fraction in models of Table 1 
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0.482 
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0.422 



extreme as M13 or NGC 2808 which we are considering here, 
so the F vs. Na anticorrelation should be a result of the pri- 
mordial abundance v ariations, which are predic ted by our 
mod els. Nevertheless, Denissenkov et al.l (j2006b) point out 
that ISmith et al.l (120051 7 data also show an average decrease 
of the fluorine and carbon abundance for increasing lumi- 
nosities of the M4 giants, and attribute the anticorrelation 
to extra mixing, explored with models similar to those we are 
presenting here. In our models with extra mixing, fluorine is 
in fact heavily depleted, especially in the models having the 
largest diffusion coefficient. Figure [5] shows the fluorine vs. 
oxygen anticorrelation and allows a prediction: in the most 
extreme giants of M13 and NGC 2808, having the lowest 
oxygen abundance which we attribute to extra mixing, the 
fluorine abundance must be more or less reduced depending 
on the fastness of extra mixing. Table 3 shows in fact that 
the diffusion coefficient is an important ingredient to achieve 
strong depletions (see the comparison between models 4 and 
5 or the comparison between models 6 and 7). 



5 THE FINAL HELIUM ABUNDANCE IN THE 
ENVELOPE 

Extra mixing aff ects the envelope helium abundance 
l| Weiss et al.l l2000l 'l. In our case, the envelope abundance 
is already a high Y=0.4, and it increases up to Y=0.48 



in the most "realistic" sequence, computed with mass loss. 
We do not explore further this problem here, but mention 
that some of our models also showed higher envelope helium 
abundances. This fact has important consequences for the 
location of these models on the following Horizontal Branch 
phases. The very hot HB stars at Toff > 31000K can not 
be explained by standard HB models, even with very low 
hydrogen envelope on top. Sweigart (1997) suggested that 
a late helium flash along the white dwarf cooling sequence 
can penetrate the small entropy barrier of the weakly burn- 
ing hydrogen-burning shell (Iben 1976), so that helium and 
carbon would come to the surface, and hotter Toff's can 
be achieved. We see in the models with extramixing that, 
maybe, the very high envelope Y is a result of deep mix- 
ing in giants which already from the beginning have a high 
helium abundance, and suggest that the existence of a hot 
blue tail in the HB may be linked to the extremely high he- 
lium abundances in the envelope, a combined result of the 
initially already high helium and of deep extra mixing. The 
evolution of these stars will be described in a subsequent 
paper (Caloi et al. 2007, in preparation). 



6 CONCLUSIONS 

The chemical anomalies in GCs include both oxygen deple- 
tion and helium enrichment. While there is no totally ac- 
cepted model for the site of nucleosynthesis which affected 
the matter from which the anomalous stars were formed, it is 
accepted that helium enrichment will go together with oxy- 
gen depletion. Helium enhancements up to Y^O.4, as seen 
in a few main sequences of GCs, may result from models, 
while the same models are easily consistent with the deple- 
tion of oxygen down to [O/Fe]~-0.4 shown by most GCs. 
The models on the contrary have serious difficulties with 
the much larger depletions which are seen in some giants 
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in M13, or may be inferred by the spectral analysis of the 
giants in NGC 2808. 

In this work we remark that the evolution along the 
red giant branch (RGB) for stars having Y~0.4 is notice- 
ably different from the evolution of normal red giants with 
Y'^ 0.24: in particular, both the milder external convec- 
tion penetration, and the higher initial helium abundance 
in the envelope conspire to produce a smaller discontinuity 
in molecular weight along the evolution (in particular, the 
standard "red giant bump", occurring when the hydrogen 
burning shell crosses the molecular weight discontinuity, is 
less prominent and occurs at a larger luminosity). 

We suggest that only the stars having very high helium 
abundance (if they are indeed present in the GC) are subject 
to deep non standard extra-mixing below the formal convec- 
tion zone, a mechanism which is not efficient in the presence 
of the molecular weight barrier. We propose that extra mix- 
ing may be efficient along the whole RGB evolution, if the 
helium content is large enough to minimize the molecular 
weight discontinuities, and attribute to it the occurrence of 
the extreme abundance anomalies which are present in a few 
GCs. 

We model extra mixing then starting from already 
anomalous chemistry, and show that it affects mainly the 
oxygen -and not the sodium- abundance. We also predict 
a very strong depletion of Fluorine, and helium envelope 
abundances as high as Y~0.48 in present models. 

These results may serve to clarify the role of primor- 
dial self-enrichment in GCs, showing that self-enrichment 
should not necessarily explain also the most drastic chemical 
compositions found in some clusters. 
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Figure A2. HR diagram. From right to left arc shown the 
0.8Mq (Y=0.24), 0.7Mq (Y=0.35) and O.64M0 (Y=0.40) 
for Z=2xl0-^. The three dots on the tracks show the location of 
the red giant bump. 



Figure A3. We show the luminosity versus time evolution of 
the three tracks (Y increasing from left to right). The times are 
arbitrarily normalized. Red giant bump of the Y=0.4 model is 
much less prominent than the others. 
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APPENDIX A: MODELLING THE EVOLUTION 
OF HELIUM RICH MODELS ALONG THE RGB 

We describe here the evolution of three masses having 
different helium content: Y=0.24 for M^O.SMq, Y=0.35 
for M=O.7M0and Y=0.4 for Y=O.4M0, and a metallic- 
ity in mass fraction Z=0.002. The masses are chosen as 
representative of masses presently evolving in GC (age of 
~ 10 — 12Gyr). We compute standard models, in which mix- 
ing is purely due to convection, and no kind of overshooting 
is allowed below the external convective region. We show in 
Fig lAll the hydrogen abundance along the RGB evolution of 
The corresponding HR diagram is shown in Fig |A2l 

Fig lAll shows that the deepening of convection in the 



Y=0.24 model soon produces a molecular weight discontinu- 
ity in the stellar envelope. The discontinuity reaches a max- 
imum at a mass fraction Mr/ Mq~ 0.27, due to the maxi- 
mum extension of convection at the so called "first dredge 
up". Also a model above the RG bump, when the burn- 
ing shell has crossed the molecular weight discontinuity, is 
shown. We see from the central and right figure that the 
discontinuity is less prominent for the Y=0.35 and Y=0.40 
models, and especially that it is not at all present dur- 
ing the first phases of evolution along the RGB. The lo- 
cation of the discontinuity is at about the same mass point 
(Mr/ Mq~ 0.27) for the three cases. We show in FiglMlthe 
HR diagram of the three tracks and the location at which the 
shell crosses the discontinuity. For the Y=0.35 and Y=0.4 
cases the effect on the luminosity evolution is barely seen in 
the tracks, and does not produce any appreciable effect on 
the luminosity function, as we can understand by looking at 
the time evolution of the luminosity, shown in Fig |A3l As 
a byproduct, notice that this result shows that we can not 
infer the existence of groups of stars very much enriched in 
initial helium content from the effect on the RGB bump. 

Coming back to the paradigm of the extra-mixing in 
the RGB evolution, we hypothesize that the extreme helium 
rich models may not be able to avoid extra-mixing along the 
RGB, and also that extra mixing may not necessarily begin 
above the bump luminosity, as there is no or a very small /i 
discontinuity which can prevent it. 



12 F. D'Antona and P. Ventura 




0,2 0,4 o,e 



Figure Al. From loft to right we siiow tlie evolution of the hydrogen profile along the three evolutions of Figure Al. The chemical 
discontinuity is very evident in the Y=0.24 case (left), much smaller in the Y=0.35 case (center) and very small in the Y=0.40 case 
(right). 



